Numerical evaluation of the Casimir interaction between cylinders 
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We numerically evaluate the Casimir interaction energy for configurations involving two perfectly 
conducting eccentric cylinders and a cylinder in front of a plane. We consider in detail several special 
cases. For quasi-concentric cylinders, we analyze the convergence of a perturbative evaluation based 
on sparse matrices. For concentric cylinders, we obtain analytically the corrections to the proximity 
force approximation up to second order, and we present an improved numerical procedure to evaluate 
the interaction energy at very small distances. Finally, we consider the configuration of a cylinder in 
front of a plane. We first show numerically that, in the appropriate limit, the Casimir energy for this 
configuration can be obtained from that of two eccentric cylinders. Then we compute the interaction 
energy at small distances, and compare the numerical results with the analytic predictions for the 
first order corrections to the proximity force approximation. 

PACS numbers: 12.20.Ds; 03.70.+k; ll.lO.-z 



I. INTRODUCTION 

It has been 60 years since Casimir [l| found a pro- 
found explanation for the retarded van der Waals in- 
teraction as a manifestation of the zero-point energy of 
the quantum electromagnetic field. For many years the 
Casimir effect was little more than a theoretical curios- 
ity. But interest in the phenomenon has blossomed in 
recent years. Experimental physicists have realized that 
the Casimir force affects the workings of micromachined 
devices, while advances in instrumentation have enabled 
the force to be measured with ever-greater accuracy. On 
theoretical grounds, considerable progress has also been 
achieved by studying the dependence of the Casimir force 
with the geometry of the conducting surfaces . 

Up to now, most experiments aiming at a measurement 
of the Casimir force have been performed with parallel 
plates 0], or with a sphere in front of a plane 01 ■ The 
parallel plates configuration has a stronger signal, but the 
main experimental difficulty is to achieve parallelism be- 
tween the plates. This problem is of course not present in 
the case of a sphere in front of a plane, but its drawback 
is that the force is several orders of magnitude smaller. 
The problem of the theoretical evaluation of the elec- 
tromagnetic force for this configuration has been solved 
recently @- 

The configuration of two eccentric cylinders have both 
experimental and theoretical interest [1, Although 
parallelism is as difficult as for the plane-plane configu- 
ration, the fact that the concentric configuration is an un- 
stable equilibrium position opens the possibility of mea- 
suring the derivative of the force using null experiments 
(for example, one could consider experimental configura- 
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tions in which a metallic wire is placed inside a larger hol- 
low cylinder). The Casimir interaction energy between 
two eccentric cylindrical shells has been computed in , 
and was initially reported in Q . Therein, it was used the 
mode summation technique combined with the argument 
theorem in order to write the Casimir energy as a con- 
tour integral in the complex plane, to end with an exact 
formula in which the vacuum energy is written in terms 
of the determinant of an infinite matrix. As a bonus, 
it has also been shown that the matrix elements in the 
general formula for two eccentric cylinders reproduce, as 
a limiting case of relevance, those of the Casimir energy 
for the cylinder-plane configuration. The latter geometry 
is also of experimental interest: being intermediate be- 
tween the sphere-plane and the plane-plane geometries, 
it can shed some light on the longstanding controversy 
about thermal corrections to the Casimir force. Keeping 
the two plates parallel has proved very difficult, while the 
sphere and plate configuration avoids this problem, the 
force is not extensive. In the case of the cylinder-plane 
configuration, it is easier to hold the cylinder parallel and 
the force results extensive in its length. There is an on- 
going experiment to measure the Casimir force for this 
configuration [loj . 

The aim of this paper is to provide a precise numer- 
ical evaluation of the analytical results obtained in 
The numerical evaluations will allow as to test differ- 
ent approximations, as the proximity force approxima- 
tion (PFA) for close surfaces, and the "around the diag- 
onal" approximation for quasi concentric cylinders. We 
will also show numerically that the energy of the eccentric 
cylinders configuration reproduces that of the cylinder- 
plane configuration, in the limit of very large eccentricity, 
when the radius of the external cylinder is also very large. 
Finally, we will present a detailed numerical evaluation 
for the vacuum energy for the cylinder-plane configura- 
tion, providing numerical support for the analytic pre- 
dictions of the first order corrections to PFA [If[ . 

The paper is organized as follows. In Section [TT] we 
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FIG. 1: Geometrical configuration for the eccentric cylinders. 
Two perfectly conducting cylinders of radii a < b, length L ^ 
a, b, and eccentricity e interact via the Casimir force. 



will discuss the general procedure used for the numerical 
evaluation of the Casimir energy in every case considered. 
In Section Hill we will evaluate the exact formula for the 
interaction energy between eccentric cylinders. The com- 
plexity of the numerical evaluation increases as the radii 
of the cylinders get closer, and we provide details of the 
size of the matrices needed to assure convergence of the 
numerical results. In Section lTVI we analyze the particular 
case of quasi-concentric cylinders. We test the validity of 
the approximation developed in Q , based on tridiagonal 
matrices, that we extend to the next order by consider- 
ing pentadiagonal matrices. In Section|V]we consider the 
particular case of concentric cylinders. We will obtain a 
new analytic result in the small distance limit, that in- 
cludes the corrections to PFA up to the second order. We 
will also present an improved numerical method to eval- 
uate the interaction energy at small distances. Finally, 
in section IVIl we will numerically show that the inter- 
action energy for the cylinder-plane configuration can be 
derived in the appropriate limit from the eccentric cylin- 
ders configuration, a result that was anticipated analyti- 
cally for the matrix elements in Q. In addition, we will 
evaluate numerically the cylinder-plane Casimir energy 
as the minimum distance between the surfaces is much 
smaller than the radius of the cylinder. We will be able 
to show numerically that the energy is well reproduced 
in this limit by the PFA, and to compute the first order 
correction to PFA for both TM and TE modes. In the 
first case (TM), the fits of the numerical data reproduce 
with high precision the analytic prediction On the 
other hand, the fits for TE modes are close to the ana- 
lytic results or not, depending on the assumption about 
the higher order corrections. 



II. NUMERICAL APPROACH 

The evaluation of the Casimir interaction energy be- 
tween two eccentric cylindrical shells (FiglT]) has been 
initially performed using PFA in Ref. [6|. However, it is 
possible to go beyond the PFA and find an exact formula 
for the interaction energy [1, Q . This can be done using 
a mode by mode summation technique combined with 



the argument theorem. By starting with the expression 
of the Casimir energy as E = {h/2) J2pi^p ~ has 
been shown in Q that the Casimir interaction energy for 
two eccentric cylinders can be written as 
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where M™{P) = det[5„p - A™] and MTE(/3) = 
det[(5„p — A™]. Here /3 is a dimensionless integration 
variable and n,p are arbitrary integers. Roughly speak- 
ing, the function M that determines the Casimir energy 
through Eq.lU]) is such that its zeros give the eigenfre- 
quencies of the geometric configuration. More precisely, 
it is the ratio of the function associated to the actual 
geometric configuration and the one associated to a con- 
figuration in which the conducting surfaces are very far 
away from each other P 
are defined as 
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and similarly, for the TE modes, we have 
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where a = &/a is the radio between the outer and inner 
cylinder's radii and e is the eccentricity (see Fig.l). /„ 
and Kn denote the modified Bessel functions. 

In order to calculate the Casimir interaction energy, 
one needs to perform a numerical evaluation of the de- 
terminants in Eq.([T]), followed by a numerical integration 
in the variable p. We find that as a approaches small val- 
ues, larger matrices are needed for ensuring convergence. 
Likewise, as a 1, the contribution to the integral is 
significative for a bigger integration range (bigger values 
of P contribute). That turns the problem into a real chal- 
lenge from the numerical point of view. We numerically 
compute the Casimir interaction energy using a Fortran 
program. Once the M matrix elements for each config- 
uration considered is defined, we use a standard routine 
to calculate its eigenvalues and determinant. Finally, we 
perform a standard integration over all values of p. The 
parameters used by the program are: the dimension of 
the M matrix (N,N), the number of addends m corre- 
sponding to each element of the M matrix, the integra- 
tion limit (/3max) and the precision desired. The difficulty 
in running the programme lays in the compromise taken 
between all the parameters chosen. 

In the following, we will evaluate the Casimir inter- 
action energy for eccentric, quasi concentric, concentric 
cylinders, and also for the particular limit of a cylinder 
in front of a plane. 
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III. ECCENTRIC CYLINDERS 

111 this section, we present the numerical results for 
the Casimir interaction energy for two eccentric cylinders 
given by Eqs.® and 
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FIG. 3: Exact Casimir interaction energy difference AE be- 
tween tire eccentric and concentric configurations as a func- 
tion of 5 = e/a for different values of a = b/a. Energies are 
measured in units of L/A-ko? . The maximum &t 5 — shows 
the instability of the concentric equilibrium position. 



FIG. 2: Exact Casimir interaction energy difference |Ai5| be- 
tween the eccentric and concentric configurations as a func- 
tion of a = b/a for different values of 5 = e/a. Here 
AE = Ei2 — Ei2. Energies are measured in units of L/Ana^ . 

In Figsll]and[31 we reproduce the exact Casimir inter- 
action energy difference A£' — E12 — between the 
eccentric and concentric configurations. In Fig|2]we plot 
the interaction energy difFerence \^E\ as a function of 
a for different values of the eccentricity 5 = e/a. These 
numerical results interpolate between the PFA and the 
asymptotic behavior for large a Q • In FiglH we show 
the Casimir interaction energy difFerence as a function of 
5 for various values of a. Again, it is evident that the 
equilibrium position [5 ~ 0) is unstable. 

In Refs.[l, Q similar plots were performed using an al- 
gebraic evaluation of the trace of the matrix M, and a 
numerical integration, both using Mathematica. Due to 
this procedure, it was possible to evaluate the vacuum 
energy only for relatively large values of the parameter 
a, in order to reach convergence. With the numerical 
method we are presenting here, we are able to include 
smaller values of a, closer to the PFA region, where pre- 
vious numerical calculations could not reach. Both in 
FiglH and Fig[3] we include runs for a > 1.75. In or- 
der to achieve so, we have used matrices of dimension 
(21,21) and 501 addends in the sums of Eqs. ^ and ^ 
to assure convergence (variation smaller than 10^"*). For 
values a > 5, smaller M matrices (5,5) can be used to 
obtain theplots with equal precision (indeed, as it was 
shown in [8|, for a 00 the energy is dominated by the 
00-element of the matrix M). The size of M for the runs 
was set by the smaller values of a that needed bigger M 
matrix to obtain the same accuracy. 

The convergence of the numerical results depends both 
on the values of a and S. For example, for a = 3 and S = 



0.01 matrices (5,5) are enough, while for 6 ~ I matrices 
(9,9) are needed. In the case of a = 1.25, it is necessary 
to use matrices of dimensions (55,55) and (101,101) when 
S is 0.01 and 0.1, respectively. 



IV. QUASI-CONCENTRIC CYLINDERS 

In this section, we consider the situation in which the 
eccentricity of the configuration is much smaller than the 
radius of the inner cylinder {6 1). For a small non- 
vanishing eccentricity, the behaviour of the Bessel func- 
tions in Eqs.® and ^ is Im-n{PS) ~ (/SJ)"""^ for small 
arguments. This suggests that the main contribution 
should be the one coming from the diagonal elements, 
and that one only needs to use matrix elements near the 
diagonal. We will test this idea through a numerical com- 
parison between the Casimir interaction energy for differ- 
ent approximations and the exact energy derived in the 
previous section. 

First order approximation. To begin with, we 
will only consider the matrix elements proportional to 
Io{PS),Ii{P6) and IiiPS) as we are assuming smaU ec- 
centricity S = e/a 1. In this particular case, the M 
matrix become tridiagonal and the e-dependent part of 
the Casimir energy will be quadratic in the eccentricity. 
We will describe in detail the case of the Dirichlet (TM) 
modes; the treatment of Neumann (TE) modes is simi- 
lar. As was already mentioned in to order 0{S^), the 
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non-vanishing elements of the matrix A™ are 
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Second order approximation. In this case, we will con- 
sider that the main contribution to the Casimir inter- 
action energy comes from the terms that contain up to 
0{5'^), extending the previous approximation to the next 
non trivial order. Then, the matrix A has additional non 
diagonal contributions, i.e. it is a pentadiagonal matrix 
with non-vanishing elements given by 
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In the following, we will numerically compare the 
Casimir interaction energy for the quasi concentric con- 
figuration computed using the first and second order ap- 
proximations described above, with the one obtained us- 
ing the exact formula given in Eqs.(l2]) and ([3]). 

In FigsHl El sin^d [S] we present the Casimir interaction 
energy diference AE = E12 — between the eccentric 
and concentric configurations as a function of 5 = e/a 
for different values of a. Therein, we can see the differ- 
ent curves corresponding to the exact energy difference 
(solid line) and the approximated ones obtained to first 
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FIG. 4: Exact Casimir interaction energy difference AE be- 
tween tlie eccentric and concentric configurations as a func- 
tion of (5 = e/a for different values of a. The solid line is the 
exact evaluation of the Casimir interaction energy, while the 
dashed line with triangles is the second order approximation 
and the dashed line with dots is the first order one. The an- 
alytic curve (dashed line with crosses) is the result of using 
Eq.(|4]). Energies are meeisured in units of L/A-rra^. 



(dashed line with dots) and second order (dashed line 
with triangles) in the eccentricity. In all cases, matrix of 
dimension (21,21) [l^l and 501 addends in the sums have 
been used to assure convergence. 

We are also comparing the numerical results with the 
"analytic" result for quasi-concentric cylinders obtained 
Ref. 3. Therein, it was shown that the determinant of 
the tridiagonal matrix can be explicitly evaluated up to 
quadratic order in S, and therefore it was possible to write 
the interaction energy as a series 
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(5) 



The contribution of the TE modes to the interaction en- 
ergy i?™ has a similar expression, replacing the Bessel 
functions by their derivatives with respect to the argu- 
ment in the equations above. The numerical evaluation 
of this formulae in the plot is presented as the "analytic" 
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FIG. 5: Exact Casimir interaction energy difference AE be- 
tween the eccentric and concentric configurations as a func- 
tion of 5 = e/a for different values of a. The solid line is 
the exact evaluation of the Casimir interaction energy, while 
the dashed line with triangles is the second order approxi- 
mation and the dashed line with dots is the first order one. 
The analytic curve is the result of using Eq.@. Energies are 
measured in units of L/Atvo^. 



0.005 



LU 

< 




FIG. 6: Exact Casimir interaction energy difference AE be- 
tween the eccentric and concentric configurations as a func- 
tion of (5 = e/a for different values of a (a smaller range of 5). 
The solid line is the exact evaluation of the Casimir interac- 
tion energy, while the dashed line with triangles is the second 
order approximation and the dashed line with dots is the first 
order one. The analytic curve (dashed line with crosses) is 
the result of using Eq.(|4]). Energies are measured in units of 
L/47ra^ 



In FigHlwe can see that the difference between several 
approximations increases with (5, for a given value of a. 
For example, for a = 5 and 6 ~ 0.5, the analytic ap- 
proximation differs from the exact result in 2%, and the 
first order approximation in 1%. For a larger value of 
(5, as (5 = 0.9, the differences are 9% and 4% for the an- 
alytic and first order approximations, respectively. The 
results for the second order approximation coincide with 
the exact result within the accuracy imposed. However, 
as the value of a becomes smaller, the differences become 
more visible. For a = 3 the accuracy of the different ap- 
proximations is, for 6 = 0.5, 12% (analytic), 5% (first 
order), and 0.6% (second order). On the other hand, for 
S = 0.9, we have 40% (analytic), 16% (first order), and 
4% (second order). 

In Fig|5]we present results for smaller values of a, and 
in Figini we show a zoom-plot in order to appreciate bet- 
ter the differences. 

In all cases considered, we can observe the expected 
hierarchy between the different approximations: while 
the second order approximation remains very similar to 
the exact result (within a 5% error) for a > 2 and S < 0.5, 
the difference between the first order approximation and 
the exact result for a = 2 and S = 0.5 is approximately 
20%. 

Taking into account that the different approximations 
are derived under the assumption (5^1, the validity of 
the approximate results is, in general, better than ex- 
pected. Moreover, the results of this section show that 
the combination of analytic and numerical results allow a 
much more efficient numerical evaluation of the Casimir 
energy. In the particular case considered here (quasi con- 
centric cylinders), from the numerical point of view it is 



much more convenient to consider sparse matrices con- 
centrated on the diagonal, than large matrices in which 
all elements are non vanishing. 



V. CONCENTRIC CYLINDERS 

In this section we will derive an analytic result for 
the vacuum energy in the concentric cylinders configura- 
tion valid for small distances, beyond PFA, and we will 
present an improved numerical method to evaluate the 
interaction energy at small distances for the particular 
case of two concentric cylinders. 

The exact formula for eccentric cylinders coincides, of 
course, with the known result for the Casimir energy for 
concentric cylinders (e = 0). Indeed, as /„-m(0) ~ (5„m, 
in this particular case the matrices A™'"^^ become di- 
agonal and the exact formula reduces to [1, [l^ : 
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The first factor corresponds to Dirichlet (TM) modes and 
the second one to Neumann (TE) modes. The concentric- 
cylinders configuration is interesting from a theoretical 
point of view, since it can be used to test analytic and 
numerical methods. It also has potential implications for 
the physics of nanotubes 0, [l^ . 



6 



The proximity limit a — 1 <IC 1 has already been an- 
alyzed for the concentric case [ll]. In order to compute 
the Casimir energy in this limit, it was necessary to use 
the uniform expansion of the Bessel functions and to per- 
form a summation over all values of n. As expected, the 
result is equal to the one obtained via the proximity ap- 
proximation, namely 
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pTE,cc _ pTM.cc _ 1 pEM.cc 
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(8) 

and both TE and TM modes contribute with the same 
weight to the total energy. 



A. Beyond proximity approximation: the next to 
next to leading order 

We will now compute analytic corrections to the PFA 
given in Eq.®. Due to the simplicity of this configura- 
tion, we will be able to obtain not only the next to leading 
order, but also the next to next to leading contribution. 
In order to do that, we need the uniform expansions of 
the Bessel functions. We have 
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where 



(10) 



the functions u{t) and v{t). They read 
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Replacing Eq. (fT2|) into Eq.®, and expanding the log- 
arithm as a series, it is possible to compute explicitly the 
remaining integrals in /?. After a long calculation, the 
Casimir energy can be written as 



« ^ — ^/i + (l + i)(a-l) 

12 360a2(a-l)3\ ^ M 4^^ ^ 

1 1 1 A / x9 



(15) 



In the expression above, the first term inside the paren- 
thesis corresponds to the proximity approximation con- 
tribution given in Eq.Q, while the second and third 
terms are the first and second order corrections, respec- 
tively. It is worth noticing that the sub-leading term 
coincides with the result obtained by means of the semi- 
classical approximation p^ . It is also important to re- 
mark that both TM and TE modes contribute with the 
same weight to the energy in the leading and the next 
to leading orders. However, this is not the case in the 
quadratic term. There is a factor I/tt^ coming from the 
TM mode, and a factor I/tt^ -1-1/10 corresponding to the 
TE one. 



B. Improving the convergence of the numerical 
evaluation 
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With these expansions at hand, we can evaluate the 
matrix M both for the TE and TM modes. The expres- 
sion in Eq. ([7]) can be approximated by 

M"" « (l - e-2"A„(y)^(„_ j _ g-2„A„(y)^^^^ ^ 

(12) 

where 



with h{y) = v/r+y2. In Eq.p^ we have defined coeffi- 
cients A{n^ y) and Bin, y) in terms of the expansions of 



Numerical calculations of the Casimir energy for a very 
close to one are very difficult since big number of terms 
have to be considered in the sums, and therefore conver- 
gence problems arise, mainly produced by underflows and 
overflows in the evaluation of Bessel functions of large or- 
ders. 

In order to perform a numerical evaluation of the 
Casimir energy in the proximity region, we will describe 
a subtraction method, in which we have used the value 
of the energy in the PFA to improve the numerics [l^. 

In the case we are concerned here, we can add and 
subtract the interaction energy for concentric cylinders 
computed using the leading uniform asymptotic expan- 
sion for the Bessel functions, up to first order in a — 1: 
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(16) 



We denote by E the interaction energy obtained by in- 
serting these expansions into Eq. ([5]), which can be com- 
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FIG. 7: Ratio between the exact Casimir energy for concentric 
cylinders E12 and the Casimir energy estimated using the PFA 
up to the next to leading order Epp^, as a function of the 
parameter a. This is done for two different methods: the 
numerical (of slow convergence) and the numerical improved 
(subtraction method). 



puted analytically 

E = - 
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2fe(a-l) 



1) 



(17) 



and contains the leading order of the Casimir energy. 
Now we write 



E- 



12 



(£;- -E) + E. 



(18) 



The difference contained in the brackets in Eq. (1181) . has 
a faster convergence than the original sum and, therefore, 
can be easily calculated numerically. 

In Fig[7] we present both Casimir energy of the con- 
centric cylinders for the direct numerical calculation (of 
slow convergence) and the alternative method mentioned 
above. In this figure we plot the ratio E12/ Epf^ where 



TT^NTL 
^PFA 



360a2(a- 1) 



1 



-(a-r 
2 



(19) 



As can be seen, with this subtraction method it is possi- 
ble to compute the exact energy for values of a much 
closer to 1, while the accuracy of the direct calcula- 
tion is worse for a < 1.02 . Moreover, the numeri- 
cal results confirm the analytic result given in Eg. ljlSp . 
We have fit the ratio between the Casimir interaction 
energy with the next to leading correction of Eg. p^ . 



-^12/^p/a = a + b{a~ ly , obtaining a = 1.00, and 
b = 0.29. The analytical results, expected from Eqs.([T5]) 
and are a ~ 1 and b ~ 0.3026, which means that we 
are checking the next to next leading correction with an 
error smaller than 1.5%. 



A similar method could in principle be applied to the 
eccentric cylinders or the cylinder-plane configurations, 
although in these cases the main difficulty is the ana- 
lytic evaluation of the approximate energy that has to be 
added and subtracted. 



VI. A CYLINDER IN FRONT OF A PLANE 





FIG. 8: Cylinder-plane configuration. A perfectly conducting 
cylinder of radius a is in front of a perfectly conducting plane 
at a distance d. 

In this section, we will study the cylinder-plane con- 
figuration (FiglH]). The Casimir energy for this config- 
uration was first evaluated in the PFA in Ref. Q. The 
exact formula has been derived in Refs.[lll, EBIj and has 
the same structure than Eq. ([1]) , but with the matrix ele- 
ments ^™ and A™ replaced by the corresponding ones 
for this geometry, that we will denote by A™'^^ and 

/iTE,CP 

71, p 

As suggested by simple geometric arguments, the ec- 
centric cylinders formula reproduces the cylinder-plane 
matrix elements in the limit 6, e — s- 00, keeping H — b — e 
fixed. Indeed, using the uniform expansion of the Bessel 
functions, and as explained in Ref.[^, the matrix ele- 
ments A™ and A™, become respectively, 



n,p 



Km 



Kn+p{2(3H/a) = X 



_ /lTE,CP 



(20) 



and 



A. 



TM 



K^+,{2fiH/a)^A™'^^. (21) 



These expressions coincide with the result for the TE and 
TM modes in the cylinder plane configuration [ll|, [l^ . 

In the following we will, firstly, compare the exact 
Casimir interaction energy for the cylinder-plane config- 
uration with that of the two eccentric cylinders, in the 
limit that the latter reproduces the former configuration. 
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In the end of the section, we will numerically evaluate 
the cylinder plane Casimir for small distances, in order 
to discuss the leading correction to the PFA. 



B. Numerical evaluation of the cylinder-plane 
Casimir energy beyond the Proximity Force 
Approximation 



a=300 
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A. Comparison between eccentric cylinders and 
cylinder-plane configurations 

We will now show explicitly that the numerical evalu- 
ation of the vacuum energy for eccentric cylinders, based 
on Eqs. (|l|2p and ([3]), reproduce the exact results for the 
cylinder-plane configuration, described by Eq.([T|) with 
the matrix elements given by Eqs. ((20|) and (|2ip . 



1.08 
1.07 
1.06 
1.05 

Q_ 

% 1.04 

LU 

1.03 
1.02 
1.01 
1 

1 2 3 4 5 6 7 8 

d/a 

FIG. 9: Comparison between the exact Casimir interaction 
energy for eccentric cylinders and the cylinder-plane config- 
uration as a function of d = H — a for different values of 



In Figl9]we present the ratio of both energies as a func- 
tion of d for different values of a. These runs were done 
by the use of matrices of dimension (81, 81) and 901 ad- 
dends in the sums of Eqs.([2l) and ([3]). The need of big 
matrices is set by the cylinder-plane configuration pro- 
gram, while the number of addends is set by the eccen- 
tric cylinders geometry. From the numerical results we 
see that, as expected, the vacuum energy of the eccentric 
cylinders tends to the vacuum energy of the cylinder in 
front of a plane for large values of h and e, when H and a 
are fixed. The coincidence is of course better for smaller 
values of d, the minimum distance between surfaces. 

In all cases, we can sec that the exact Casimir inter- 
action for eccentric cylinders is bigger than the cylinder 
plane energy. This result is expected from the PFA, since 
the conducting surfaces are closer to each other in the 
case of the two eccentric cylinders than in the cylinder- 
plane configuration, for a given minimum distance be- 
tween surfaces. 



In this section we present a detailed computation of 
the vacuum energy for the cylinder-plane configuration. 
In Fig llOl we present the Casimir interaction energy for 
the cylinder-plane configuration obtained by the use of 
our Fortran program. For the runs, wc used a matrix of 
dimension (101,101) to reach the proximity area {d 0). 
It must be mentioned that for smaller values of d, we 
need to increase the dimension of the A matrix and the 
integration range of (3 in Eq.(IT]). This fact becomes our 
major limitation to reach yet smaller values of d. 




FIG. 10: Leading term for the Casimir interaction energy of 
the cylinder-plane configuration. The leading term behaves 
proportional to -0.0228 d"^/^. A simple fit f{x) = -yx" of the 
numerical data for d/a <1 gives 7 = —0.021 and i, = —2.53. 



Wc now discuss in more detail the limit d <^ a. This 
problem has been considered from an analytical point 
of view in Ref . [ll]| . Using the uniform expansions for 
the Besscl functions appearing in the matrix elements 
^np'^^ and A™'^^, and after complex calculations, it 
can be shown that, in the proximity limit: 



^CP 



a 3C(4) / 

rf^ 32V2 V 



0.1944- 



27r V d^ 32 V2 V a 



(22) 



(23) 



where we have written separately the contributions of 
TM and TE modes. 
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FIG. 11; Numerical result for the TM modes for the cylinder- 
plane configuration, and the corresponding fits presented in 
Table in A simple linear fit f{x) = a + bx of the numerical 
data in the interval 0.04 < d/a < 0.07 gives a = 0.9999 and 
b = 0.1900. The theoretical values are a = 1 and b = 0.1944. 



d/a 


/i(x-) = l + bx 


.f2{x) — l + b*x + c*x^ 


f'i{^) = l + b*x + c*x^* log(a;) 


[0.04 : 0.15] 


b = 0.1864 


b^ 0.1922, -0.0601 


6 = 0.1961, c = 0.0438 


[0.04 : 0.20] 


b = 0.1849 


6 = 0.1923, c= -0.0613 


6 = 0.1983, c = 0.0540 


[0.04 : 0.25] 


b = 0.1829 


fe = 0.1922, c= -0.0601 


b = 0.2003, c = 0.0634 


[0.04 : 0.30] 


b = 0.1811 


6 = 0.1920, c= -0.0586 


b = 0.2022,c = 0.0716 


[0.04 : 0.35] 


b = 0.1794 


6 = 0.1918, c= -0.0572 


6 = 0.2045, c = 0.0810 


[0.04 : 0.40] 


b = 0.1771 


6 = 0.1914, c = -0.0549 


6 = 0.2076, c = 0.0935 



TABLE I: Different fits for the numerical results of Fig. 1111 We fix fi{0) = 1 since the numerical data agree this value with 
high precision. 



d/a 


Mx) = l+bx 


.f2{x) — l + b*x + c*x'^ 


f[i{x) = l + b*x + c*x^* log(a;) 


[0.04 : 0.15] 


b = -0.8301 


fe = -0.9704, c= 1.4499 


b = -1.0711, c= -1.0852 


[0.04 : 0.20] 


b = -0.8013 


fe = -0.9509, c= 1.2326 


b = -1.0772, c= -1.1141 


[0.04 : 0.25] 


b = -0.7683 


6 = -0.9349, c= 1.0794 


b = -1.0890, c= -1.1674 


[0.04 : 0.30] 


b = -0.7399 


6 = -0.9222, c = 0.9772 


fe = -1.1037, c= -1.2306 


[0.04 : 0.35] 


b = -0.7158 


6 = -0.9091, c = 0.8879 


6 = -1.1232, c= -1.3115 


[0.04 : 0.40] 


b = -0.6851 


6 = -0.8943, c = 0.7999 


b = -1.1534, c= -1.4360 



TABLE IL Different fits for the numerical resuhs of Fig. \T2\ 
high precision. 

We will discuss the first order corrections to PFA for 
TM and TE modes separately. In Fig llli we show our nu- 
merical results for the TM modes. The fit of the numer- 
ical results depends of course on the interval chosen for 
d/a. There is an obvious compromise: on the one hand, 
as already mentioned, we cannot consider very small val- 
ues for d/a because of numerical limitations. On the 
other hand, the expansion in powers of d/a are expected 
to be valid only for d/a <C 1. In any case, as can be seen 
from Table I, the different fits for the numerical results 
are stable, and confirm both the PFA to leading and next 
to leading orders. Indeed, the results are fully compati- 



We fix /i(0) — 1 since the numerical data agree this value with 

ble with the analytic results given in Eg. ip^ . considering 
both linear and quadratic fits of the numerical results. 
Moreover, a simple linear fit in a smaller range of d/a 
gives a = 0.9999 and b = 0.1900 and already reproduces 
the analytical results [ill with high accuracy (see also 
numerical findings in [l7[). 

In Fig |121 we show our results for the Neumann modes, 
and we include in Table II different fits of the numeri- 
cal data. In this case, the value obtained for the linear 
correction to PFA depends strongly on the assumption 
about the next non trivial correction. This is not sur- 
prising: as we cannot consider extremely small values for 
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FIG. 12: Numerical result for the TE modes for the cylinder- 
plane configuration, and the corresponding fits presented in 
Table UT] A simple linear fit f{x) — a + bx of the numerical 
data in the interval 0.04 < d/a < 0.07 gives a = 0.9940 and 
b = —0.7808. The theoretical values are o = 1 and b = 
-1.1565. 
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FIG. 13: A numerical fit of the results for the TE modes 
including cubic corrections /(x) = 1+bx+cx^ \ogx+dx^ . The 
coefficients are b = -1.0478, c = -0.9485, and d = 0.6708. 

d/a, the non linear corrections may have a non negligible 
contribution in the intervals chosen for the fits. For ex- 
ample, a simple linear fit gives a = 0.994 and b = —0.7808 
which does not coincide with the result in Eq. (j^H]) . How- 
ever, based on the discussion about the slower conver- 
gence of the Neumann corrections presented in Ref . pH . 
we have allowed the possibility of non linear corrections 
proportional to (d / a)'^ hi(d / a) in our fits. Remarkably, 
when this non linear corrections are taken into account, 
the coefficient of the linear correction gets closer to the 
analytic prediction in Ea. (|23p . that we reproduce with an 
error less than 7%. Note that, as can be seen in FigfTTl 



this is not the case for TM modes, since the best fit of 
the numerical data contains a quadratic term without a 
logarithm. In Fig. 13 we show a fit of the numerical data 
for TM modes that includes a cubic correction (d/a)^. 
With this additional term, the fit reproduces the numer- 
ical data up to d/a = 0.5. 

To summarize, the fits of the numerical data clearly 
confirm the analytic prediction for the TM modes, 
and suggest that the next non trivial correction for 
the TE modes is not quadratic but proportional to 
{d/af\n{d/a). 



VII. FINAL REMARKS 

We have numerically evaluated the Casimir interaction 
energy for the two eccentric cylinders configuration and 
for the cylinder plane geometry, extending in several di- 
rections the numerical results presented in Ref.[^. For 
quasi concentric cylinders, we have shown that the ap- 
proximation based on tridiagonal matrices derived in [8| 
is in good agreement with the numerical values. We also 
extended this approximation to the case of pentadiagonal 
matrices. Our results show that, for small eccentricities, 
it is far more efficient to consider the contribution of the 
matrix elements near the diagonal, than a "tour de force" 
numerical calculation based on the exact formula. 

For concentric cylinders, we have obtained analytically 
the quadratic corrections to the PEA. As far as we know, 
this is the first explicit non linear correction to PFA ex- 
isting in the literature. We have also shown that the PFA 
can be used as a useful tool in order to improve the nu- 
merical evaluation at very small distances, and we have 
used this improvement in order to check numerically the 
non linear correction to PFA. 

Finally, we have analyzed in detail some numerical re- 
sults for the cylinder-plane geometry. On the one hand, 
we have shown that the Casimir energy for this configura- 
tion can be obtained from that of the two eccentric cylin- 
ders. Although this coincidence has been anticipated for 
the matrix elements in Ref. the numerical data show 
that the result is also valid for the energy. On the other 
hand, we have computed the TE and TM contributions 
to the energy for small distances, and compared the fits 
of the numerical results with existing analytic predictions 
for the linear corrections to the PFA. 
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